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ANNEALED INVARIANCE PRINCIPLE FOR RANDOM WALKS ON RANDOM 
GRAPHS GENERATED BY POINT PROCESSES IN R d 

ARNAUD ROUSSELLE* 


Abstract. We consider simple random walks on random graphs embedded in R“ and generated by point 
processes such as Delaunay triangulations, Gabriel graphs and the creek-crossing graphs. Under suitable 
assumptions on the point process, we show an annealed invariance principle for these random walks. These 
results hold for a large variety of point processes including Poisson point processes, Matern cluster and 
Matern hardcore processes which have respectively clustering and repulsiveness properties. The proof 
relies on the use the process of the environment seen from the particle. It allows to reconstruct the original 
process as an additive functional of a Markovian process under the annealed measure. 


Key words: Random walk in random environment; Delaunay triangulation; Gabriel graph; point 
process; annealed invariance principle; environment seen from the particle; electrical network. 

AMS 2010 Subject Classification : Primary: 60K37, 60D05; secondary: 60G55; 05C81, 60F17. 


1. Introduction and results 


1.1. Introduction. Random walks in random environments are natural objects for studying flows, molec¬ 
ular diffusions, heat conduction or other problems from statistical mechanics. Many papers concerning 
random walks on random perturbations of the lattice Z rf or random subgraphs of Z rf were published dur¬ 
ing the last forty years. Particular examples include random walks on supercritical percolation clusters 
and the so-called random conductance model (see Bisl 1 ] and references therein). An annealed invariance 
principle for such random walks was obtained in [DMFGW85 DMFGW89] while its quenched counterpart 
is proved in SS04 BB07 BP07, MP07, Mat08| . Few papers were written in the case where the random 
environment is generated by a point process in M d . In the last decade, techniques were developed to 
study random walks on complete graphs generated by point processes with jump probability which is a 
decreasing function of the distance between points possibly modified by the presence of an energy mark 


at each point. The main results on this model include an annealed invariance principle FSBS06,FM08 


isoperimetric inequalities CF09 , recurrence or transience results CFG09 , a quenched invariance princi¬ 
ple and heat kernel estimates CFP13|. In |Roul5b|, we answer the question of recurrence or transience of 
random walks on random graphs constructed according to the geometry of a random (infinite and locally 
finite) set of points such that Delaunay triangulations or Gabriel graphs. In the present paper, we are 
interested in stating annealed invariance principles for simple random walks on random graphs generated 
by point processes. More precisely, given a realization £ of a point process in M rf , a connected graph 
G(£) = (£, Eq(q), in which each vertex has finite degree, is constructed. While its vertex set is £ itself, 
its edge set E G ^ is obtained using deterministic rules based on the geometry of £. For instance, we will 
consider Delaunay triangulations, Gabriel graphs or creek-crossing graphs generated by point processes. 
Their constructions will be described below. We then consider continuous-time random walks on G(£) 
and investigate their annealed (or averaged) scaling limits. 
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A. ROUSSELLE 


In the sequel, we will write A f (resp. Mo) for the collection of all infinite locally finite subsets of R d 
(resp. infinite locally finite subsets of containing the origin). We denote by V (resp. Vo) the law 
of a stationary, ergodic and simple point process (resp. its Palm version). We recall that Vq(Mo) = 1 
(see |SW08 §3.3]). We denote by £ (resp. So) the expectation with respect to V (resp. Vo). We assume 
that P-almost all £ is aperiodic, i.e., with ^-probability 1: 

£ + Fz£, Vx G R d , 

where t x is the translation r x £ = Sy-x- This implies that the same holds for "Po-almost all £°. 

For almost any realization of £°, we define a variable speed random walk (VSRW) on G(£°) as follows. 
Let D(M + ;£°) be the space of cadlag £°-valued functions on M + endowed with the Skorohod topology, 
and xf = X^°(t), t > 0 be the coordinate map from D( R + ;£°) to £°. For x G £°, we write P x 

for the probability measure on Cl^o := -D(M_|_;£ U ) under which the coordinate process (Xf )t>o is the 

£° ~ 
continuous-time Markov process starting from AT = x with infinitesimal generator given by: 


^°f(y) ■= <£*(/(*) ~ y e 


( 1 ) 




£° £° 
where Cy, z is equal to the indicator function of the event ‘{y,z\ G Eqi^o}. We denote by E x the 

expectation with respect to P x . Note that P x can be seen as a probability measure on := Z?(M + ; 
concentrated on Sl^o. 


*) 


£° 

The behavior of X% under P x can be described as follows: starting from the point x, the particle 
waits for an exponential time of parameter deg G ^oj(x) and then chooses, uniformly at random, one of its 
neighbors in G(£°) and moves to it. This procedure is then iterated with independent times (of parameter 
the degree of the current location in G(£ 0 )). The same arguments as in FSBS06, Appendix A] show that 


these random walks are well-defined (i.e. that no explosion occurs) as soon as £o[deg G (go)(0)] < oo. 
Let P be the so-called annealed semi-direct product measure on A/"o x Cl: 


P[P(£°,A«°)] = f Pf[F(e, AT°)]P 0 (d£°). 

J No 


£° 

Note that (X f s )t >o is no longer Markovian under P. Using technics similar to those from 


DMFGW89 


Section 4] and |FSBS06j, we will use the point of view of the particle and [DMFGW89, Theorem 2.2] to 


obtain annealed invariance principles for random walks on specific examples of random geometric graphs 
with suitable underlying point processes. Let us describe the precise models. 

1.2. Conditions on the point process. In all what follows we assume that V is the distribution of a 
simple, stationary point process £ satisfying £[#(£ n [0, l] d ) 8 ] < oo. We also suppose that the second 
factorial moment measure of the point process £: 

A (2) (Ai x A 2 ) := £ [#{(xi,x 2 ) : x\ G £n A u x 2 G £n A 2 \ {xi}}], A U A 2 G B(R d ), 

is absolutely continuous, with bounded density, with respect to the 2d-Lebesgue measure. We assume that 
V is aperiodic, isotropic (i.e. its distribution is invariant under rotations around the origin) and almost 
surely in general position (i.e. there are no d + 1 points (resp. d + 2 points) in a (d — l)-dimensional 
affine subspace (resp. in a sphere), see |Zes08j). Using e.g. [CFP13, Lemma B.2], one can see that the 
Palm version Vq of V also satisfies these properties. It is moreover assumed that £ has a finite range of 
dependence k, i.e., for any disjoint Borel sets A,BcM d with d (A, B) := inf{||x — y\\ : x G A, y G B} > k, 
£ n A and £nB are independent. In particular, the point process is ergodic. 

Finally, the following hypothesis on the void probabilities (V) and the deviation probabilities (D) are 
needed: 
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(V) there exists a positive constant c\ such that for L large enough: 


V 


#(^n[—L/2,L/2] d ) =0 


< e 


—c\L d . 


(D) there exist positive constants 02,03 such that for L large enough: 


V 


# (en[-L/2,L/2] d ) >c 2 V 


< e 


-c 3 L d 


Assumptions (V) and (D) imply the following analogue estimates for the Palm version of the point 
process (see Lemma HU 0 and ([ 8 ])): 

(V’) there exists a positive constant C 4 such that for L large enough: 


Vo 


#(e°n(x + [-L/2,L/2] d )) =0 


< e 


-c A L d 


\/x G 


(O’) there exist positive constants cs,C 6 such that for L large enough: 


Vo 


#(e°n(x+[-L/2,L/2f)) 


>c 5 L° 


< e 


-c 6 L d 


\/x G 


Let us note that, due to the stationarity of the point process and (D), we have: 

lim N & V\#(£ PI [— N, N] d ) > ciN d 1 = 0 for some 07 > 0. 

A—>oc L v J 

These assumptions are in particular satisfied if £ is: 

• a homogeneous Poisson point process (PPP), 

• a Matern cluster process (MCP), 

• a Matern hardcore process I or II (MHP I/II). 

We refer to Appendix[B]for an overview on each of these point processes. Note that these processes have 
different interaction properties: for PPPs there is no interaction between points, MCPs exhibit clustering 
effects whereas points in MHPs repel each other. 

1.3. The graph structures. Given a realization £ of the point process satisfying the previous assump¬ 
tions, we write 

Vor^(x) := {x G R d : ||x — x|| < \\x — y||,Vy G £} 

for the Voronoi cell of x G £; x is called the nucleus or the seed of the cell. The Voronoi di agra m of £ 
is the collection of the Voronoi cells. It tessellates into convex polyhedra. See M0l94,CallO| for an 
overview on these tessellations. 

The graphs considered in the sequel are: 

DT(£): the Delaunay triangulation of £. It is the dual graph of its Voronoi tiling. It has £ as 
vertex set and there is an edge between x and y in DT(£) if Vor^(x) and Vor^(y) share a (d — 1)- 
dimensional face. Another useful characterization of DT(£) is the following: a simplex A is a cell 
of DT(£) iff its circumscribed sphere has no point of £ in its interior. Note that this triangulation 
is well defined since £ is assumed to be in general position. 

These graphs are widely used in many fields such as astrophysics jRBFNOl , cellular biology 
[Pou04], ecology |Roq97] and telecommunications |BB09 


G n (£), n > 2: the creek-crossing graphs generated by £. These are sub-graphs of DT(£). Each 
G n (£) has £ as vertex set and two vertices x, y G £ are connected by an edge in G n (£) iff there do 
not exist an integer k < n and vertices uq = x,..., uj. = y G £ such that || U{ — Uj+i|| < ||x — y|| 
for all i G {0,..., k — 1}. Let us recall that G 2 (£) is also known as the relative neighborhood 
graph generated by £ and that the intersection of the G n (£) (in the graph sense) is the Euclidean 
minimum spanning forest of £ (see [HNS13, HNGS15]). 
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Gab(£): the Gabriel graph of £. Its vertex set is £ and there is an edge between x, y £ £ if the ball 
of diameter [x, y] contains no point of £ in its interior. Note that Gab(£) is a subgraph of DT(£) 
and contains the Euclidean minimum spanning forest of £ and the relative neighborhood graph. 
It has for example applications in geography, routing strategies, biology or tumor growth analysis 
(see [BBD02, GS69,MS80|). 

1.4. Main results. The aim of this paper is to prove the following theorem. 


Theorem 1. Assume that V satisfies the conditions described in Subsection 1.2. 

Let (Xf ) t > o be the random walk on DT(£°), Gab(£°) or G n (£°), n > 2 with generator 
Then, the rescaled process Yfi° ,£ := {eXl_ 2j )t >o converges weakly in Vo-probability to a nondegenerate 
Brownian motion W a ~ ld with: 


= ton ]So [£f [0*f • e i) 2 ]] > °> 


( 2 ) 


where e\ is the first vector of the standard basis of M d . In other words, for any bounded continuous 
function F on the path space D([0, , oo), M d ) endowed with Skorohod topology 

A 0 




(p° 


E 


F W a ~ ld 


in Vo-probability. 


This is the case if the underlying point process is a PPP, a MCP or a MHP I/II. 


In the companion paper |Roul5a 


we prove a quenched invariance principle in the particular case of the 
simple random walk on the Delaunay triangulation generated by suitable point processes (see also |FGG12 


for the 2-dimensional case). The method used in |Roul5a| relies notably on heat-kernel estimates and 
demands a precise control of regularity (isoperimetry, volume growth, ...) of the random environment. 
The main interest of the present paper is that it does not require such precise estimates. Consequently, 
it applies to a wider class of point processes and graphs. 

As mentioned above, we closely follow the method of | FSBS06 and we want to apply the results of De 
Masi, Ferrari, Goldstein and Wick. A key idea of DMFGW89 , which also appears in the prior works by 
Kozlov Koz85 and by Kipnis and Varadhan |KV86|, is to introduce an auxiliary Markov process: the 


- - ... . £° 

environment seen from the particle. Actually, thanks to the aperiodicity of the medium, (Xf )t> o can be 

reconstructed as a functional of the environment process. 

We will prove the following preliminary results stated in the general context of the introduction which 
imply Theorem [l For x,y £ £°, we write x ~ y if {x, y} £ Fg^o). This definition depends on £° which 


is omitted from the notation. 

Proposition 2. Let V be a stationary aperiodic ergodic simple point process in M. d . Assume that 
max x ^o ||x|| and deg G ^o)(0) admit polynomial moments of any order under Vo- 

Then, the rescaled process Y}°' £ := {eX^_ 2t ) t >o converges weakly in Vo-probability to a Brownian motion 
W D with covariance matrix D characterized by equation (19). 

Remark 3. (1) If the underlying point process is isotropic, then D = <r 2 I^ where a 2 is given by Q). 

(2) In the applications presented here rnax x ^o ||x|| and deg G ^o)(0) actually admit exponential moments 

Mil¬ 


as shown in r Roul5a 


Note that the diffusion matrix D of W D can be degenerate. We will show that, under additional 
assumptions, it is not. To this end, we will need to discretize the space and to consider boxes of side 
K > 1 : 


= B z 


K 


= Kz + 


K Kid 


-— , z £ Z d . 

2 2 . 1 ’ 


We will also define so-called ‘ good boxes' in the same spirit as in |Roul5b| which allow us to deduce ‘nice 
properties’ of G(£) from ‘nice properties’ of a suitable percolation process on 7L d . 
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Proposition 4. Let V be the law of a point process satisfying the assumptions of Proposition^ Assume 
in addition that £[#(£ n [0, l] d ) 8 ] < oo and that there exists c 7 > 0 such that: 

lim jV 8 p[#(£n [-N,N] d ) < c 7 N d ] = 0. (3) 


If one can associate to V-a.a. f a subset of the boxes {B z , z G Z rf } ; called (£-)good boxes, such that: 

(1) in each good box B z , one can choose a reference vertex v z € B z n 

(2) there exists L such that to each pair of neighboring good boxes B Zl and B Z2 , one can associate 
a path in G(£), (v Zl ,... ,v Z2 ), between the respective reference vertices v Zl and v Z2 of these boxes 
satisfying: 

(a) (v Zl ,..., v z ,,) C B Zl U B Z2 , 

(b) the length of ( v Zl ,... ,v Z2 ) in the graph distance is bounded by L, 

(3) if f is distributed according to V, the process X = { X z , z G Z d } := {1 b z is £-good, z £ Z d } 
stochastically dominates an independent Bernoulli site percolation process on h d with parameter 
P>pf e ( Z 2 ), 

then the limiting Brownian motion W D in Proposition [^] has a nondegenerate covariance matrix. 


1.5. Outline of the paper. Sections [2] and [3] are respectively devoted to the proof of Propositions [2] and 
|4| Theorem [l] is derived from Propositions [2] and [4] in Section [4| More precisely, arguments which allow to 
deduce polynomial moments for deg PT( -g O) (0) and max x ^ 0 in dt(£°) ll x ll under Vo from the assumptions on 
the point process are given in Subsection [4T] In Section 4.2, we explain how to construct the ‘good boxes’ 
of Proposition [4] in the case of creek-crossing graphs. The nondegeneracy of the limiting Brownian motion 
then follows in the other cases by Rayleigh monotonicity principle. Appendix [B] contains the verification 
of the fact that PPPs, MCPs and MHPs satisfy the assumptions of Theorem [l] 


2. Proof of Proposition I 


We show how to use IDMFGW89, Theorem 2.2] in our framework. This result was used in the context 


of random walks in random conductance models (see DMFGW85, DMFGW89]) and of random walks 


on complete graphs generated by point processes with jump rates depending on energy marks and point 
distances (see |FSBS06 ). It allows us to obtain the convergence of the rescaled process to a Brownian 
motion in Proposition [2] and gives an expression for the diffusion matrix of the limiting Brownian motion. 
The nondegeneracy of this matrix and the proof of Proposition [4] are postponed to Section [3] For this 
section, we assume that the assumptions of Proposition [2] are satisfied. 


2.1. The point of view of the particle. Let us denote by (T x ),j. gR[ i the group of translations in 
which acts naturally on N as follows: t x £ = by-x■ Note that, if £° 6 No and x G £°, then t x £ 0 G No- 
The theory developed in |KV86,DMFGW89| makes extensive use of the process of the environment seen 
from the particle (r f o£°)t>o- Whereas the original random walk is defined on the probability space 


(0^o,Pq ), this auxiliary process is defined on Z° := D(M + ;A/"o) and we write Pto for its distribution. 
A generic element of is denoted by = (£^ )t>o- The environment process can be seen from an 
‘annealed’ point of view, that is under the distribution P := f P^o[-]'Po(d^°). In this setting, under 
the assumption of aperiodicity, this is the continuous-time Markov process with initial measure Vo and 
transition probabilities given by: 


p[£WI£ = f] =p*[tf = ?] : = 


p o [ X t = x ] if C = r x ^,xef Vs t > 0 
0 otherwise 


















6 


A. ROUSSELLE 


Let us explain how the original walk can be reconstructed from the process £°. For t > 0, let Xt : 
be the random variable defined by: 

M?) ■■= E ( 4 ) 

s e [0, t] : 

A s (|°)/0 

where A s (£°) = y if there exists y such that = T y ^°_ and is 0 otherwise. Note that A s is well defined 

for a.a. realization of the environment process since is almost surely aperiodic. The law of (Xt(£°))t>o 

f° ~ 

under P^o is then equal to Pq , the distribution of the original random walk starting at 0. In particular, 

the diffusion matrix can be calculated using the following equality: 

£ 0 [E^°[(xf -a) 2 ]} =E[(X t -a) 2 }. 

The family := Xt — X s }t >iS >o is an antisymmetric additive covariant family of random variables 

in the sense of DMFGW89], and X t has cadlag paths. Hence, to apply |DMFGW89, Theorem 2.2], we 
only need to check that: 

(1) the environment process (Ct)t> 0 is reversible and ergodic w.r.t. P; 

(2) the random variables t > s > 0 are in L 1 (H°, P); 

(3) the mean forward velocity: 

exists as a strong limit in L 1 (Mo 1 'Po)', 

(4) the martingale M t := X t — <p(£®)ds is square integrable under P. 

The rest of this section is devoted to the verification of these hypotheses. 

2.2. Underlying discrete-time Markov chains and description of the dynamics. It is convenient 
to decompose the dynamics in terms of the jump chain and the jump times as described in this subsection. 
Given £° E Mo, let us write := (£°) N and denote by (x| )„. e N a generic element of For x E £°, let 

pf be the distribution on of the discrete-time simple nearest neighbor random walk on G(£°). We 

denote by Ex the corresponding expectation. On another probability space (0, Q), let Tn,x, x E n E N 
be independent exponential random variables with mean 1/ deg G ^o)(x), and define on (fl^o x 0, pf <g) Q) 
the random variables: 

n— 1 


K 




:=0;P« U :=^T«° 0 ,n>l; 


»=o 


i, Xf 


nf (t) := n where R < t < P^ +1 . 

These random variables can be thought respectively as the instant of nth jump and the number of jumps 
up to time t for the continuous-time random walk with environment starting at the origin. The same 
arguments as in |FSBS06, Appendix A] show that almost surely no explosion occurs, i.e. 

Pj° <g> Q \ lim P|° = oo] = 1 . 

L nfoo J 

n to •—'tO 

In particular, for a.a. £ , n* (t) is well defined for all t > 0 and (X ? 0 ) t >o is a jump process which 

n% (t) 

satisfies the same dynamics (under pf <8) Q) as (xf )t >o (under pf ): starting from the point x, the 
particle waits for an exponential time of parameter deg G ^o)(x), chooses uniformly one of its neighbors in 
G(£°) and moves to it. 

It is useful to introduce the probability measure Qo on Mo given by: 

n (A^\ - deg G(«°)(°) v 

2o(d? ’ - fo[deg G(s o )( 0)] P » (d « >' 
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and P := f Pq Qo(d£°), the annealed law of the (original) discrete-time random walk. In the same way 
as in Subsection 2.3, one can obtain reversibility and ergodicity of the discrete-time process under P. Let 


E denote the expectation under P. 

One can also consider the discrete-time version of the environment process (Cn)neN := (F~ c o£ 0 ) ne N and 
write P^o for its distribution on the path space := ATq. When is aperiodic the same identification 
as in Subsection 2.1 is possible. It is easy to see that the processes (deg G ^o)(x| )) ne pj on (H^o,Pq ) and 
(deg G ^ O s(0)) ree N on (H°, P^o) have the same law. One can also verify that the process (l^) n eN dehned by 
Yq = 0 and: 

n—1 

Y n :=^2 A (^ 0 ,^ +1 ), 

i =0 

with A(g ) ,g ) +1 ) defined in the same way as in (jdj), has distribution Pq on fl^o. If there does not exist 
such that = r x £, we set A(£,£') := 0. 


Remark 5. By slight modifications of the arguments given in Subsection 2.3, one can obtain reversibility 
and ergodicity of (£°)neN under P := /P^oQo(d^ 0 )- This can be used to derive an annealed invariance 
principle for the discrete-time random walk thanks to ]DMFGW89 , Theorem 2.1]. 


2.3. Reversibility and ergodicity. 

2.3.1. Reversibility. We want to show that the environment process (fit )t>o is reversible with respect 
to P and we follow the argu men ts of the proof of JFSBS06 , Proposition 2]. Using the description of the 
dynamics given in Subsection 2.2, one can check that, for any n E N and any £o = £°, £i, • • •, = C Q/ E A/q: 


nf (t) = ra, C% = 6, • • •, C% = 

ri^ Jrt n 


a 0 ' 


nf ( t ) = n, o =£ n - C% = Co 


Ft' 


Ft: 


and deduce that: 


P e o [$ = C°'] = P^o' [4° = C°], *>0. 


Hence, for any positive measurable functions f,g and any t > 0, thanks to Neveu’s exchange formula 
(see |SW08| Theorem 3.4.5]), we obtain: 


E[/(Co°MC?)] = / E/^ 0 )^ r ^ 0 ) p ?°fe 0 = ^ 0 ]^o(dC 0 ) 


'No 


l No 


'No 


xe£° 


E Pr^O [e? = C°]^o(dC°) 




E /(lA(C°) P ? ofe° = T,C 0 ]^o(dC 0 ) 


ce£° 


= E[/(C°MCo°)]- 

The reversibility of (C?)t>o then follows. 


2.3.2. Ergodicity. We follow the proof of |FSBS06, Proposition 2]. Note that it suffices to ascertain that 
for any measurable A C Mq satisfying P^ofC^ E A] = 1a(C°) for a.a. C°) we have Vo (A) = 0 or 1. Given 
such an A, there is a subset B C A such that Tq[A \ B] = 0 and P^ofC? E A] = 1 for all C° E B. For 
x E C°, C° E B , one then has P^ofC? = t x C° E A\ = P^ofC? = T x f°] > 0 by the positivity of the jump rates. 
Thanks to FSBS06, Lemma 1 (iii)], V$(A) = 0 or 1 which completes the proof. 
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2.4. Integrability of X[ S)t ] = Xt — X s , t > s > 0. This is an immediate consequence of: 

Lemma 6 . Let 7 G N* = {1, 2,...} and to > 0. There exists a constant 0 < c 8 = c 8 (j, to) < 00 such 
that: 

£0 [Eq [ll^f II 7 ]] <c 8 Vt, for all t< t 0 - (5) 


To obtain this lemma the reader should redo the calculations of the proof of FSBS06, Proposition 
1 ] with p = q = 2 using that deg G ^o)( 0 ) and max x ^o INI admit polynomial moments and make the 
dependence on t apparent. 

2.5. Existence of the mean forward velocity. We want to verify that the mean forward velocity: 

exists as a strong limit in L ] . Actually, we prove the following stronger result which is also useful for the 
square integrability of t 

Lemma 7. Let ip : A/"o 


square integrability of the martingale Mt in Subsection 2.6 

R d defined by: 


- E 


= 2-, x - 


( 6 ) 


Then, 


Proof. As in Subsection 
to 




2.2 


we denote by nf (t) the number of jumps of the environment process starting 


at £ u up to time t. First, we use that Xtl^o Q = Xq = 0 and (a + b) 2 < 2(a 2 + b 2 ) to write: 


1 2 XVo 


|E e o [X t \ -Mf) || 2 Po(de°) 


? -M« 0 ) + E c»Ai nf( , te2 ]|| 2 n(dO 

L ll E ? [ x ‘h ?* VK»)|| 2 A(df») + !/ 11%7,i nf(lfe2 ]H 2 A(df). (7) 


2 

~t 2 -Wo 


We need to show that the two summands in ([7]) go to 0 with t. We first deal with the first term in the 
r.h.s. of (7]). 

Note that, with the notation of Subsection |2.2| for x G £°: 

P € o = r x e,nf{t) = 1] < pf <g> Q[xf = x, t |° 0 < t\ 

= P? [Xf = x](l - e- des W°) {0)t ) 

-Ss77M degGK ° ,(0)f = ll ~ oi ’ 

where we used that 1 — e~ u < u. Thus, l x ~ot ~ P^o [f,t = (t) = l] is nonnegative. 

It follows from the Cauchy-Schwarz inequality that: 


£b[||E£o[X t l n4 o =1 ] -f<^°)|| 2 ] = Sol’ll E ^(P?ofe° = r x f,nf(t) = l] - l^o t) 

xG£° 

<£o[S 1 (e)S 2 (f°)}, 

Si(f) ■■= E (l^ot-P e ofe 0 = r^°,nf(f) = l]), 


2i 


with 


(8) 
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and 

‘ S ' 2 (?°) := E IMI^xMrf-PfO [£° = r x t?,n?(t) = 1]). 

Let us observe that 

Si(£°) = deg G(? 0 ) (0)t - P ? o [nf (f) = l] 

= deg G (^o)(0)t - P^o [nf (t) > l] + Pjo [nf (t) > 2] 

= deg G(? o)(0)t - 1 + e _deg G«°) (0 )* + P ?0 [nf (t) > 2] 

< deg G(? 0 ) ( 0 ) 2 t 2 + P ? o [nf (t) > 2], 
where we used that for n > 0, 0 < — 1 + e~ u + u < u 2 . Moreover: 

P^° [ni (t) > 2] < Pq < 8 > Q [Tq j0 < t, T^ o < t] 

= (1 - e - de s G « 0 )(°)t) ^(1- e~ deg G(^? 0 )(°)f Pco g = r^°] 
xe£° 

< t 2 deg G(e o) ( 0 )% [ deg G( ^ i) ( 0 )], 

thus, with ([£]) and ( 10 ): 

Si (£°) < t 2 ( deg GK o) ( 0) 2 + deg G(? o) ( 0 )E 5 o [ deg G( £ } (0)]). 

On the other hand, one has: 

= E ||®|| 2 (lx~oi - P$o g = T x £°,nf(t) = l]) < *E IM| 2 - 

:r~0 


Note that: 


£0 


d egc«o)(0) 2 E 


rr~0 


< £0 


deg G(? o)( 0 ) f 


‘So 


max \\x\ 


Thanks to the definition of P in Subsection 2. 2l we can write 


Tn 


deg G(? o ) (0)E ? o[deg G( ^ ) (0)] E ||x|| 

tc~0 

= £0 [ deg G( ^ 0) (0)] E deg G( ^ } (0) ||x| 




< ^0 [ deg G(€ o) (0)] E [ deg G(fi) (0) 2 ] 5 E [ ( E 


tc~0 


2-i 


Due to the stationarity of the discrete-time environment process under P, one has: 


So 


deg G (£ 0 )( 0 )E ? o [deg G(i r } ( 0 )] E ||x|| 

<£ a [ deg G( £ 0) (0)] 5 E deg G(s o } (0) 


‘So 


EH 2 de gG«o)(°) 




VI 

deg GK o } ( 0) 3 

~ 2 S 0 ' 

to 

VI 

deg GK o } ( 0) 3 

~ 2 £o~ 


1 

4 ^ 




Thanks to ([ 8 ]), (11), (12), (13) and ( fl4| ), the first summand of ([7]) goes to 0 as t —> 0. 
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Let us move to the second term in the r.h.s. of Q. Applying the Holder inequality with exponents 4 
and 4/3 and the Cauchy-Schwarz inequality, we obtain: 

121 <£:o[E ? o[||A f || 4 ]^P ? o[nf(t)>2] 


£o 


l E ?° e° 


(t)> 2 J 


ll 1 


< £o 


E ? o[||Ai 


on 


P € o[nf(t)>2] : 


One can obtain as in (10) the alternative upper bound: 

Pfo[raf (t) > 2] < fE^o [ deg G( |'o ) (0)]. 


(15) 


(16) 


Hence, (10) multiplied (16) to the square, the Jensen’s inequality and the stationarity of the discrete-time 
environment process under P = J Pgo<2o(d£°) imply that: 


£c 


P 5 ° [nf(t)> 2] ; 


< t 4 £ 0 

< t 4 £o 

.4 


d e g G (£°) (0)E^o [ deg G( ^o ) (0)] 3 

deg G (^o)(0)E € o (deg G( |'o ) (0)) 
= i 4 £o[deg G( £ O) (0)]E (deg^o^O)) 3 
= t 4 £o \ (deg G( ^ 0) (0))' 


(17) 


Finally, with (J5|, ( |15| ) and ( j 1 7[ ) ? wg obtain that for t smalh 

^£o[||E € o[Xa 4 0 (t) J|| 2 ] <cgti, 

which tends to 0 with t. □ 

2.6. Square integrability of the martingale M t = X t — J ^ y>(£°)ds. Due to LemmaJfjJ X t E L 2 (H°,P), 
so it is enough to see that: 

/%(4°)d S GL 2 (S°,P). 

Jo 

Using the Cauchy-Schwarz inequality and the stationarity of the environment process under P, one has: 


E 


r W°)ds" 21 


Jo 


<tE\ [ \\v(0\\ 2 d S ]=t [ E[||^ s 0 )|| 2 ]ds = t 2 E[||^ 0 0 )|| 2 ] 
L Jo J Jo 

= * 2 £o[|M| 2 ] < t 2 fo[deg G ^o)(0) 4 ] 2f 0 [(max||x||) 4 ]2 < oo. 


2.7. Infinitesimal square displacement as a L 2 -limit. By following precisely the same proof as in 
Lemma [7] with higher order moments estimates, one obtains: 

Lemma 8. For E A/o, let if) be the d x d symmetric matrix characterized by: 

(a • i/’(^°)a) := ^(a • x) 2 , a E M rf . (18) 




Then, (a ■ Ei 2 (A(o,Po) and: 


Tc v \ 21 l2 (Xo,Vo) / ,t\\ 

-E .{(a-Xt) J -——> (a • ip(-)a). 

Since its assumptions were verified in this section, Theorem 2.2 of [DMFGW89] implies Proposition 
[2] It gives moreover the following characterization of the diffusion matrix D of the limiting Brownian 
motion: 

roc 

(a ■ Da) = £0 [(a • ^a)] — 2 / (ip ■ a, e tC p ■ a)-p 0 , a E R rf , (19) 

Jo 

where £ is the generator of the environment process and ip , ip are defined in ([6]) and (18) respectively. 
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3. Proof of Proposition [4j the diffusion coefficient is nondegenerate 

From now on, we assume that the hypotheses of Proposition [4] are satisfied and we prove it implies the 
non-degeneracy of the limiting Brownian motion. 

3.1. Variational formula. In order to show that the diffusion matrix of the limiting Brownian motion 
in Proposition [2] is nondegenerate, it is useful to obtain a variational characterization of this matrix. The 
use of variational formulas to approximate diffusion coefficients is rather classical. It appears in particular 
in papers based on the works of Kozlov Koz85 or De Masi, Ferrari, Goldstein and Wick |DMFGW89 


In principle, it allows to derive the nondegeneracy of the diffusion coefficient and to obtain numerical 
estimates for it. See FSBS06,Bisll and references therein. 

Lemma 9. The diffusion matrix of the limiting Brownian motion in Proposition [1] is given by: 

(a ■ Da) = inf f £ 4° x ((a • x) + V*/(£ 0 )) V( 0 (d£°), a G M d , 

/eL°°(V 0 ,P 0 ) 


( 20 ) 


1 


where V x is the discrete gradient V x f(£,°) := /(t x £°) — /(£°). 

Proof. Minor changes to the proof of |FSBS06, Proposition 3] show that the generator C of the environ¬ 
ment process is non-positive and self-adjoint with core Vo). Moreover, for any / G L°°(J\fo,Vo), 

one has: 

Af(£°) = E c C(/(^ 0 )-/(^°))’ 

x££° 

where Cg x equals the indicator function of the event ‘{0,x} G E G ^od. We refer for example to 
Chapter X] for basic tools to prove this. 

Hence, following the same strategy as in FSBS06 and invoking equations (45) and (46) of this paper, 
one obtains: 

(ip ■ a,e tC ip ■ a)-p 0 dt = sup (2(<p ■ a, f) Vo - (/, (-£)/)?>„). (21) 


RS75 


f 

J o 


sup 

/eL°°(VoW>) 

On the other hand, Neveu’s exchange formula gives: 


E 4 0 j(f) 2 we 


'Vo 




'Vo 


'Vo 


E c T 0 fj(r x e) 2 V o(d£°) 


a;e?° 


E 4Lf(TxX°) 2 we 


xe£° 


thus, 


•'■Vo ^cfO 


x€£° 

CO 


E 4*(/(^) 2 + /(^°) 2 - 2 f(T X e)f(n°))v 0 (de 

x££° 

X <€(w/R°)) 2 n« 0 ). 

In the same way, Neveu’s exchange formula gives: 

[ E <€> • x )f(^)V o(d?°) = - [ E c oU a • x )f( T ^°)E o(d£°) 

J A/n J A/n 


2 JjVo 


'Vo 


(22) 




i£U 
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which implies that: 


(<P ' a. f)v 0 = / E c £( a ' z)/(£°)Po(d£ 0 ) 

= - 9 / o ( d ^ 0 

w ■'JVo _,n 


( 23 ) 


Finally, with equations (19), (21)-{23), and the definition of ip, one has: 


(a • Da) = £0 [(a • ipa)\ - 2 sup (2 (ip • a, /)^ 0 - (/, (-£)f) Vo ) 

feL°°(Af 0 ,r 0 ) v 7 

= f r ^ ( f ° [(° ' ^°)] “ ' °> + 2 (/> (-Qf)p 0 ) 

f£L™(jVo,Po) V / 

in f . / E C S,{(a-x) 2 + 2(a-x)V a ,/(e°) + (V a; /(e 0 )) 2 }Po(de 0 ) 

°(Afn.Vn) J \r„ z 1 J 


f£L™(jVa,r a ) J ATq 




which gives the required result. □ 

For sake of simplicity, we assume from now on that the underlying point process is isotropic. Never¬ 
theless, the arguments below can be easily adapted to the anisotropic case. 

3.2. Cut-off on the transition rates and lower bound for the diffusion matrix. In the sequel, 
we will consider the (largest) subgraph of G(£°) in which each edge has length at most r c > 0. For r c > 0 
fixed, we set: 


— 


-x,y 


= 1 


{x,y} is an edge of G(£°) and \\y— aj|<r c — Fr,y 


Thanks to the variational formula (20), we immediately obtain that: 
(a ■ Da) > 


inf / 

feL°°(Ar 0 ,Vo) Jj\f Q 


EO^ + Va/te 0 )) 2 ? o(d£°). 


(24) 


(25) 


xG£° 


As above, one can show that the operator C defined by: 

£f(?) ■= E C(/(^°) - /(£°))> / e v@), 

has core L°°(J\fo,Vo), is self-adjoint and is the generator of a Markov process (£°)t >0 on D(M + ; A/o). We 
write P ? o (resp. P) for the probability measure on D(M_|_;A/"o) associated with this Markov process with 
initial distribution 5^o (resp. Vq). The same arguments as in Subsection 
is actually equal to: 

^ poo ^ 

(i a ■ Da) := £q [(a ■ ip a)] — 2 / (<p • a, e tc Cp ■ a)-p Q dt, 


show that the r.h.s. of (25) 


(26) 


where 


£(£°) ; = E and (° • ^°) a ) : = E €> ■ x ) 2 - 

xe£° 


Hence, (25) is equivalent to the inequality D > D. 

3.3. Periodic approximants and electrical networks. The aim of this subsection is to bound from 
below D by using diffusion coefficients of random walks in appropriate periodic electrical networks. The 
method is very similar to the one in |DMFGW89,FSBS06 and relies on the theory of electrical networks 
(see |SW08 LP14| ) and homogenization (see | JK094||0wh Q31 ). Actually, we will define suitable random 
resistor networks with diffusion coefficient which can be expressed as an effective conductance and which 
is lower than D. In the sequel, the medium will be periodized along the first coordinate direction. 
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2 N 



Q points of TjY 
% points of T]\f 
+ points of i 1% 

- edges of 0y which belong to ^(0 

- edges of 0 y which do not belong to C(£) 

— ■ edges of C(0 which do not belong to 


Figure 1. The vertex set of 0^ contains the points of = £ n C^ag and Tjy. An 
edge of G(£) between points of 0.^ N is an edge of 0^ iff its Euclidean length is smaller 
than r c . There is also an edge in 0^ between each point of £ n B jy (resp. £ n B^) and 
each point of FJj (resp. Tjy). 


3.3.1. Periodized medium. For £ G AT a locally finite subset of and N > r c , we will construct two finite 
networks from G(£). To this end, let us define the following subsets of 

C 2N :=[-N,N ] d , C 2N :=]-N,N [ d , Q^:=£nC 2N , 

:= {x = (xi,... ,Xd) G z d : x\= ±N, \xj\ < N, j = 2,..., d}, 

Tljy := 0^- U Tjy U r^-, Bjy := {i £ C 2 n '■ £] — AT, — + r c ]}, 

5+ := {x G C 27V : xi G [AT - r c , AT[}, 5^ := £2^ D 

A finite (unoriented) graph = (23^r, 0^r) is then constructed; it has an edge between x and y if: 

• i,i /6 0.^ N , {x, y} is an edge of G(£) and \\x — y\\ < r c , 

• x G and y G T“ , 

• x G -E>|) + and y G Tjy. 

Another graph 0^ = (23^, 0^) is obtained from 0^ by identifying vertices: 

x- = (~N,x 2 , ■■■,x d ) and x + = (N,x 2 , ■ ■ -,x d ). 

We denote by p the projection map from Tllr onto Note that p ln e is nothing but the identity. 

The graph 0y r will represent the periodized medium in which any vertex in Tat := p(r^-)(= p(Tjy)) is 
connected to any vertex in B^ N := B^f~ U B^ + . 

We consider now a continuous-time random walk on 23^ with infinitesimal generator given by: 

£at/Of) : = c v(T>y)(/(y) - f( x )), VxG23^, 
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where 


4fi x ,y) ■= < 


if {x, y} £ <£% with x, y £ £3^ 
if {re, y} £ with i or y G Tjv . 


(27) 


#(IV) 

„ 0 if {x, y} fL 

Since C N (x,y) = c^ N (y,x), the generator C^ N is self-adjoint in L 2 (2 3^ N ,m^ N ) where m^ N is the uniform 
distribution on 93^. In general &^ N is not connected and is not ergodic. In this case, ergodic 
probability measures are uniform distributions on the connected components of and one can apply 
the results of De Masi, Ferrari, Goldstein and Wick in each connected component. We denote by (resp. 
Pj^ x , x £ 93^) the distribution on the path space := Z3(M+; 93^) of the random walk with generator 
Y and initial distribution m^ N (resp. 5 X ). In order to apply one more time 
let us introduce d\ the antisymmetric function on 93]^ such that: 

2/1 -xi if x, y £ £1% 

di{x, y) ■■= < yi + N if y £ £3^, y\ < 0 and x £ V N 
yi — N if y £ £3^, yi > 0 and x £ Tn 

Moreover, for t > 0, u = (w s ) s >o £ £1%, we define the random variable: 


DMFGW89 


Theorem 2.2], 




L i/, ) /(^ : = Y di(u a -,u a ). 

s e [o, t] 

lj _ ^4 uj s 

The last quantity is the sum of the increments of the projections along the first coordinate of the positions 
taken by the particle up to time t. It is then possible to apply DMFGW89, Theorem 2.2] to the time- 
covariant antisymmetric family := — X^^} t>s >o and obtain in particular: 

Lemma 10. For N £ N, N > r c and £ £ A I: 


1 . 


^ E At[(^AT,i’ ) 2 ] = D N := N ‘Pn)tb.% < ^', 

where tp^ N and are the scalar functions defined by: 

<Pn( x ) : = Y c N( x ^)M x ^y) and if%(x) := ^ c%(x, y)di(x, yf. 


(28) 






3.3.2. Relationship between D and the diffusion constants of periodized media. We compare D (and there¬ 
fore D) with suitable averages of the diffusion coefficients of the periodized media. 


Lemma 11. It holds that: 


and 


lim £ [m^(v4)] = £o [Vh,i 


N—>oo 


lim £[{ip%,e tc ^^ N ) i 

iV—>• OO N 


(29) 


(30) 


where (p\ (resp. is the first coordinate element of (p (resp. the first diagonal element off:) defined 

in Subsection \3.S\ 

In particular, 

"W (31) 


D > D > ( lim sup £ [D 

^ TV—>-00 
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Proof. We first derive (31) from (29) and (30). By the isotropy of Vq. D is a multiple of the identity 


matrix. Hence, using (|26|) and Fatou’s lemma: 

roo 

D > D = (£ n l^i,i| — 2 / ((p u e tL Lpi)'p Q dt)l d 

roo 


(So $i,i 

= (^ £ [’4(4)] - 2 / ^£[<4, dt )^ 

v u 

/ /*°° c \ 

> ( lim sup £ [m^( 4 )] - 21 iininf / £[( 4 > e ^ 4 ) m i 4 ) I(i 


N—>oo 


= ( lim sup £ [D ? v ] ) l d . 
' N—¥oo ' 


cO cO ° 

We now prove (29). Recall the definition of ci,y from (24) and note that c| )3/ = 0 if x G CW- 2 ?- c and 

O 

V 0 CW- Thus, for x G O.N-r c - 

^i,iijxO = ( e i ' ^{TxC) e i) = X l{o,?/} is an edge of G(t x £), ||y||<r c ^i 

y&-r x £ 

= X {vi-xif = ^%{x). 

y e *4 
{*.!/} 6 

Since a: G O.N-r c , its neighbors in 6 /v are points of £3 at and thus the last sum runs only on points of Qn- 
This allows us to write 


m 


N 


( 4 ) = 


i 


#(£ n c 2 at) + #Tat 


X] $l,i(t-*0 + X] 40) 

x£ < S^ N \Q^ N _r 




and to conclude by taking the P-expectations in this identity and by using Lemma [l 6 ] ([5]) and Lemma 
[^with h^ N = ippj. Before moving to the proof of (30), let us note that the same method shows that for 
1 < p < 4: 

lim £[n4(|</4| p )] = £o[|£i| p ] < r p c £ 0 [#(£° D B 2 (0, r c ) p ] < oo. 


N—>oo 


In order to prove (30), let us fix 0 < a < 1, set M := N — [1V Q J and define the exit times 
4 := inf {t > 0 : u t 0 C 2 N- 2 r c }, u = {ut)t> 0 6 ^ = T>(M + ; TJ^). 

Since (e tC Nip^(x) = E^ x [^(cut)], we can write: 

£ [(4> e * £]V 4) m « r ] = £ [^1,AT + 4TV + 4, n] > 

where 

A^ n '■= m N(^xeC2M^N( x )^‘N,x[^T^<t^N( UJ t)] S j, 

^3 ,N := m A r ( : *-3;eC2M4( x )4,x[ : *-r i $>t4( a; i)])- 

Then (30) will be established once we show: 

lim £[A\ n ] =0, i = 1,2 and lim £[A^] = (£i,e*4i)7V 

_/V—^OO ’ /V—Vnn ’ 


(32) 

(33) 


TV—>-oo 


( 34 ) 

















16 


A. ROUSSELLE 


For the first limit in (34), we use the Cauchy-Schwarz inequality and the stationarity of the process uj_ 
under to write: 

mi — £ m A r (,^-x&C2M | ( a ') [| ) |] ) 


< £ 

< £ 

= £ 


m 


TV 


(24 \ C 2M ) 


'£ 


tn 


N 


( 1 Pn( x ) 2 '®N,x kW 2 ]) 


+ #(£ n c 2 n \ c 2 m ) 

#^n + #(£ n c 2 n) 

#r n + #(£ n c 2 n \ c 2 m ) 
#rv + #(£nc 27V ) 


£ 


m 


•jv 


£ 


m 


■N 


{Vn( x )*) 


(^Ar(^) 4 ) 


£ 


m 


N 


^ E AT,rc [^Ar( W i) 4 ]) 


The first factor of the r.h.s. goes to 0 thanks to the dominated convergence theorem. Indeed, the quotien 
is bounded by 1 and goes a.s. to 0 since its numerator is of order IV 0 ^ 1 L-/V“J while its denominator is of 
order N d . The second factor remains bounded thanks to (32). Let us turn our attention to the second 
limit in 


Applying twice the Cauchy-Schwarz inequality and using stationarity, one has: 
l^[^ 2 ,A r ]| — £ m A r (lxeC 2 M |^At(t) I < ^A r (^) |] ) 

— ^ [ m lv l^) 2 ) ] 2 £ m N{^-xeC 2 M^ l N,x[^T^<t i £N( UJ t) 2 ]) 

- ^[ m Af(^Af( x ) 2 )] 2 < ^[ m Ar(^'A,a:[v ; Ar( a; ) 4 ])] J< ^ m N [^N — *] ) 

= ^[ 4 (^(s;) 2 )]^[rn^(^(s) 4 )]^[m^(ixec^P^[T| < t])] 4 . 

, the first two terms are uniformly bounded. The last factor goes to 0 by 


Once again, thanks to 
Lemma [l 8 l 

To deal with the last limit in (34), we need to introduce the following exit times: 

Tn,x '■= inf{t > 0 : x + X t (£°) 0 C 2 N- 2rc }, C 2 m , £° = (it)t> 0 £ S°, 


(35) 


with Xt(£°) as defined in (4) and as defined in Subsection 


2.1 


Note that, if N is large enough (i.e. [A""] > r c ), one has for x G C 2 m- 

Vn( x ) = fii(. T xO and E^Jl jt >t <P%(utj] = E r ^ [l7w,*>t£i(£?)] ■ 

Hence, 

m A (^-x£C2M l Pl 


£I a I,n\ =£ 


In the same spirit as what we did to obtain lirriv^oc A 1 ] = we S e ^ using Lemma 
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(47): 


m A \ ^-x€C2M l Pl(' T x£)'^T x £ [^-TN,x<t i P'i-{it 


so that 


lim £ 

TV—>-00 


lim £ L4| M \ = lim £ 

iV-> 00 L ,S ’ JVJ A-> 00 


= 0, 


m N (<Pi [£1 (i\ 


Finally, we apply Lemma the function 
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This function is in L 2 (A/’o,T’o) since the process is stationary w.r.t. E and deg G ^o)( 0 ) admits polynomial 
moments under TV This gives 




N—foo 


□ 


3.3.3. Relationship between D^ N and the effective conductance of suitable resistor networks. We want to es¬ 
tablish a link between the diffusion constant D^- and the effective conductance of an appropriate electrical 


network. To this end, recall definitions of Subsection 


3.3.1 


4r0, y) := c%(p(x),p(y)), {x, y} G 


— A 


and equip edges of with conductances: 


where is given by (27) and p is the projection map from 33^ onto 03^. Then, the effective conductivity 
k^ n of this electrical network is defined as the total amount of current flowing from Tjy to Tjy w hen a 
unit potential difference is imposed between these two sets. From Dirichlet’s principle (see LP14 §2.4] 
for example), we know that: 


= \ inf ' 


S ^ 4r(*>y)( v (*,i/)5 ) 2 : 9\r- =0, 5 |r + = 


= 1 


(36) 


y: (a:,2/}ee^ 


where V^g := g(y) - g{x). 

The same method as the one used in Subsection 
Lemma 10) satisfies the variational formula: 

1 


3.1 


shows that the diffusion coefficient D\ t (defined in 




, inf 

/■ Mi¬ 


ll E 


c^ r (x,y)(d 1 (x,y) + V ^)/) 2 


1 




inf 


/ : an 


Y J2 c %( x i y) (V(i,j,)(7ti + f))‘ 


(37) 


/ constant on p 1 ({x}),\/x G x£$3 N y. {x,y}££pj 


where n\ is the projection along the first coordinate. 
We claim that: 


D n ~ 


- 7 - inf 

f ■ 


f\r~ ur + — N 
11 2V Ui N 


E E_ 

seQJk v■ {x,y}e£% 


l N 


( x ,y)(V(x, y )(Ki + f)Y 


(38) 


To see this, let us fix the values of such a function / on £2^-, say = h lCl z . We will verify that the 


‘|QJ 


function minimizing (37) subject to this condition is actually constant on T^y U Tjy. Note that, when the 
values of / on £f N are fixed, we know the value of: 

and, writing C x for the value of / on p _ 1 ({x}), x E Tw, we are led to minimize in C x : 


Y | ( 2/1 +N + h(y) -c x ) 2 + Y (y 

l y eB^~ y&B^ + 


i-N + h(y)-C x f 


x£r N K y&B ' 

The summand depends on x only through C x and we must choose C x as the (unique) minimizer of this 
quantity for each x G Tjv. Thus the second infimum in (37) is achieved for functions which are constant 
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on T n U Tjy. Since / appears in the variational formula in a gradient form, one can assume that / is 
equal to N on U and (38) is proved. 

Now, let us note that the functions g : i —y M satisfying g| r - = 0 and (?| r + = 1 can be expressed as: 

1 


9 = 2 iV^ 1 + / )’ 




where / is a function on QJh such that /i r -, , r + = N and vice versa. Hence, 

\ L N UL N 

4 N 2 




#{^n) 


inf 


4r( ;r >y)( V (x,y)ff ) 2 : 9 : ^ 


9\r 


- = 0 , 


9 \rj 


= 1 


x£W% y: {x,y}££% 


Thanks to Dirichlet’s principle (36), we have just proved: 

Lemma 12. 




8 N z 






Lower bounds on are obtained in Subsection 


3.4 


3.4. Lower bound on k^ n . In order to obtain lower bounds on k^ n for large N, we use the discretization 
of the space and the ‘good boxes' introduced in Proposition [4j These ‘good boxes’ will allow us to 
link ‘nice properties’ of the random graphs we consider with ‘nice properties’ of independent Bernoulli 
processes on Z d . The rest of the proof then relies on techniques from percolation theory which also 


appear m 


FSBS06, Section 6 ]. This method is actually quite classical and is partially described in the 


2-dimensional case in Kes82, Chapter 11], 


3.4.1. Percolation estimates and LR-crossings. Let us consider an independent Bernoulli site percolation 
process Y = {Y z , z E Z rf } on 7L d with parameter p > pf te ( Z 2 ). We will use fine estimates about the 
number of left-right crossings (LR-crossings) of R 2 N, 2 (N-i) := [~N,N] x [-N + 1,N - l] d_1 . Let us 
recall that a LR-crossing with length k — 1 of Rj 2 n, 2 {n-i) f° r a realization of {Y z , z e Z d } is a sequence 
of distinct points (zi,..., z*) in i? 2 7 V, 2 (jV-i) C Z d such that |zj — Zj + i| = 1 for 1 < i < k — 1, Y Zi = 1 for 


1 < i < k, z^ = —N, z^ = N, —N < zf } < N for 1 < i < k 


and zp = zj 1 for any l > 3 and for 
4 < i < j < k. Two crossings are said to be disjoint if all the involved zjs are distinct. Note that, by 
definition, LR-crossings are contained into 2-dinrensional slices. The techniques of |[Gri89, Sections 2.6 
and 11.3] can be used, in the context of site percolation, to derive: 


(0 


Proposition 13. If p > p c (2), there exist positive constants a = a(p), b = b(p), and c = c{p) such that 
for all NSW: 

P[Y has less than bN d ~ 1 disjoint LR-crossings in R 2 N, 2 (N-i)] — ce~ aN . (39) 

Thanks to the domination assumption ([3]) in Proposition [4j there exists a coupling of the process of 
the good boxes X and the supercritical independent Benoulli procces Y on Z rf such that almost surely 
Y z < X z for all z E Z d . In other words, any good box corresponds almost surely to an open site in the 
percolation process Y. 

From now on, we fix r c := \fd + 3 K. Consider a LR-crossing '-yP erco = (zi,...,Zfc) of R. 2 N, 2 (N-\)- 

One can associate to 7 P erco a LR-crossing j & u:n Q f &^ KN from a point of F^^ to a point of r ^ N - To 
construct this LR-crossing, we first link A'zi E to the reference vertex v Z2 of B Z2 . Since v Z2 E B Z2 , 

\\Kzi — v Z2 1| < r c and {Kz\,v Z2 } is an edge of Next, for i = 2 ,..., k — 2 , we go from v Zi to v Zi+1 

using the path given in ([2]) of the definition of the good boxes in Proposition |4j These edges are edges 

of thanks to the choice of r c . Finally, we link v Zk _ 1 to Kz^ E F^ A? . We assume that 7® KN is 

self-avoiding; if not, we remove loops. Since the boxes B Zl ,..., B Zk are good, it is easy to see that: 
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(a) L( 7 ®^ jv ) < L x L( 7 perco ) where L(-) denotes the graph length of a path, 

(b) 7 ®™ c U zG 7 p-co {B z n C 2 kn), 
in particular, 


(c) if 7 perco and 72 CIC ° are two disjoint LR-crossings of R 2 n, 2 (n~i) f° r 


then 7 ®- KN and 7 ® XJV are 




disjoint LR-crossings of © 


N- 


3.4.2. Lower bound for for N-good configurations. We say that a configuration £ e AT is N-good 


if the corresponding realization of the percolation process Y has at least bN d 1 disjoint LR-crossings 
of R 2 n, 2 (n~i)- From now on, we assume that £ is an iV-good configuration and we write Cat(£) for a 
collec tion of at least bN d ~ l such disjoint crossings. By the arguments given at the end of Subsection 
there are at least bN d ~ l LR-crossings of C 2 kn hi the graph We write as before 7 ® KN for 


3.4.1 


the LR-crossing in ($^ N associated with the LR-crossing 7 P erco of R 2 n, 2 (n- 1 ) f° r Y. Since the paths are 
assumed to be self-avoiding and the 7 P erco s in C/v(£) are disjoints, it holds that: 

L( 7 ®™) < L L( 7 perco ) < L(2N + l) d . 

7 p erc °eCjv (0 7 perco e Cjv (0 

With Jensen’s inequality, one can deduce: 

1 . #(Ca (£)) 2 ^ {bN d ~^ 2 


> 


- > 


T 


7 perco gCjv (^) 

with cio = c\o(d,K). 


^ L( 7 *4») ' E 7 p..««„(5) L(7®*») ' L ( 2N + 


> c w (KN) 


d—2 


(40) 


We are now ready to bound from below n^ KN . To this end, let us introduce another electrical network 
®kn with vertex set ^ KN := U (F kn^kn}’ where f ^N’^kn are two extra points. Its edge set 


^kn consists in the collection of edges of with both endpoints in 0 ^^- and extra edges between 

. We assign a conductance 1 to each of these edges and we define the effective 


f^ iV and each point of B' 
conductance /cl~ r of ©| 


■£>± 

KN 


KN '-’kn as the total amount of current flowing from F KN to f f <N when a unit difference 

of potential is imposed. This new network can be seen as obtained from &^xn by merging together vertices 
in T 


KN and V+ N respectively. Since the same difference of potential was imposed between T KN and F^ 


KN 

in the dehnition of i^ KN , it is not difficult to see that f^ KN = L^kn- 

The disjoint LR-crossing 7 &KN of ©/ GAr exhibited above can be used to bound h^ KN from below. They 
correspond to disjoint (except at their endpoints) paths from F^^, to in the new network. By 

Rayleigh monotonicity principle, one can obtain a lower bound on f^ KN by removing each edge which does 
not appear in these crossings. Thanks to the series law each of these crossings 7 ® KN has ‘total’ resistance 


L( 7 ® ifJV ). Finally, due to parallels law and (1401), one has: 


^kn ~ k kn - 


E 

7 perco eCjv(£) 


1 


> ciq(KN) 


d—2 


L(y ®kjv) 

3.5. Conclusion by collecting bounds. Due to ( |31[ ), it remains to show that: 

lim sup 8 [-D^] > 0 . 


(41) 


TV—)>oo 


But, with Lemma 12 and (41), we obtain: 


lim sup 8 > 8 K 2 lim sup 8 


N—>oo 


TV—^oo 


N 2 




,#(2T 


liv) 


KN 


> cn lim 8 
N —>-oo 


{KN) d 
kn) 




4 is A r -good 
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where we used that: 


£ 


(. KN) 2 ^ 

—rw r K R-iv-*-^ is v-bad 
tFK^kn) 


> 0. 


Observe that thanks to (39): 


£ 


(. KN)° 


kn) 


4 is Af-bad 


= £ 


(. KN) d 

#(«na 2K „) + #(r 2A -„) 1 « is "- b » d 

< c 12 N d V[fi is iV-bad] < c 13 N d e~ aN 


N—yoo 


A 0. 


Thus, we obtain that: 


lirn sup £ ( D%] > cn lim £ 

N—toc N^oo 


(. KN) d 
# (^a'tv) . 


> 0 


thanks to Lemma m- Hence, the diffusion coefficient of the limiting Brownian motion in Proposition 
[2] is nondegenerate and Proposition [4] is proved. 


4. Derivation of Theorem Q] from Propositions [2] and [4] 

In this section, we show that the hypotheses listed in Subsection |1.2| imply the assumptions of Proposi¬ 
tions [2] and [4] in the cases of Delaunay triangulations, Gabriel graphs and creek-crossing graphs. Theorem 
[I] then immediately follows from these propositions. 


4.1. Polynomial moments for deg DT ^o)(0) and rriax x ^ 0 ; n dt(£°)(0) under Vq. In order to apply 
Proposition [2] to random walks on Delaunay triangulations, Gabriel graphs and creek-crossing graphs, we 
only need to show that the assumptions of Subsection |1.2| imply that deg G ^o)(0) and max x ^o(0) admit 
polynomial moment of any order under Vq. Since the Gabriel graph and the creek-crossing graphs are 
sub-graphs of the Delaunay triangulation, it is enough to do so in the case of Delaunay triangulations. 

Lemma 14. If the point process satisfies (V 7 ) and P’), then deg DT /£cn (0) and max x ^ 0 in dt(£°)(0 ) admit 
polynomial moments of any order under Vo- 

Proof. Let us consider a collection {C^,j € Jl} of cubes of side L that, as in Figure [2J surround Cl = 
[—(2 \y/d\ + 1 )L, (2 \y/d\ + 1 )L] d . Note that HJl = 0(L d_1 ). Let Al be the event ‘#(6*^ n £°) > 1 for 
all j G Jl and if (Cl fl £°) < C5(8[v / d] + 2 ) d L d ' where: 


C L 


— (41" V~d\ + 1 )L, (41" y/d\ + 1)L 


On Al, since each Cff contains a point of £°, points on the boundary of Cl are within distance at most 

\fdL from the nuclei of there respective Voronoi cells. Note that 0 £ can not be one of these nuclei. It 
follows that the union of these Voronoi cells separates 0 (which is in its own cell) from the exterior of the 
union of the cells intersecting Cl- Thus, on Al, we have: 

degDT(£°)(0) — f (^l H (°) < C5(8[~\/cf| +2 ) d L d . 
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Figure 2. If each (grey) cube of side L that border [— (2\Vd] + l)L,{2\\fd\ + 1 )L] d 
contains at least a point of £°, then the nuclei of the Voronoi cells intersecting the boundary 
of [—(2 1 ~\fd\ + 1 )L, (2\Vd] + 1 )L] d lie in the region delimited by the red curves. In this 
case, the union of these cells surround 0 G £° and the degree of 0 in DT(£°) is bounded by 

#([-(4\Vd] + l)L,(A\Vd] + l)L] d ne). 


Renee, for c 5 (8\Vd] +2) d L d < D < c 5 (8\Vd] + 2) d (L + l) d , (V’), (D 7 ) and the union bound imply 
that: 


Vo 


degDT(£°)(0) > V) 


< Vo deg DT ^ 0 )( 0 ) > 05(81 " Vd] + 2 ) d L c ‘ 

< Vo [.A c l } 

< Y Vo [# {cf n £°) = 0 ] 


jeJ L 


+ Vo [# (C L n £°) > c 5 (8fv^l + 2 ) d L° 

< 0{L d ~ l )e~ CiIjd + e ~ c e(s\Vd]+2) d L d = o(e~ cl4Ld ). 


Finally, for k £ N*, we have: 


£0 


(deg DT (£0)(0)^ < Dk Vo deg DT (go)(0) > D 

J D =1 

Lc 5 (8rVrfl+2) d j 

< Y^ D k V 0 deg DT ^ 0 )(0) > D 

d =1 

OO 

+ Yh 0 (L d<Kk+1 ^ l )Vo deg DT ^ 0 ^(0) > C5(8|"v / d] + 2 ) d L a 

L =1 

Lc s (8r%/rfH-2) d J 00 

< Y ° k + 0(L d(fc+1) - 1 e" Cl4id ) < 00 . 


D =1 


L =1 


Polynomial moments for max^^Q j n DT ^ 0 )( 0 ) can be proved to be finite similarly noticing that, on Al, 
max irv0 j n DT (^o)( 0 ) is generously bounded by Vd{^\\/~d\ + 1)L. □ 




























































22 


A. ROUSSELLE 


4.2. Good boxes and paths constructions for applications of Proposition [4} Due to the sub¬ 
graph relation G n (£) C Gab(£) C DT(£), it is enough to construct good boxes and paths for the creek¬ 
crossing graphs G n , n > 2 to apply Proposition [4] and obtain Theorem [lj Note that the constructions can 
be easily made by hand in the case of Delaunay triangulations as in the Chapter 2 of |Roul4 . 

We will use [HNGS15 Lemma 8 ] to verify assumptions of Proposition |4j 


Lemma 15 (Lemma 8 in |HNGS15]). Fix n > 2. Assume that £ is distributed according to a simple, 
stationary, isotropic and m—dependent point process satisfying (V) and (D). Assume in addition that its 
second factorial moment measure is absolutely continuous with respect to 2d—Lebesgue measure and its 
density is bounded from above by some positive constant. 

There exist then 0 < C 15 , ci 6 < 00 such that: 


V [- 4 g„,x,l] > 1 - C15 exp ( - L Cl6 ), 

where Ag u , x ,l is the event ‘any two vertices of G n u,v £ x + are connected by a G n -path 

contained in x + [-L, L] d \ 


In the sequel, we assume that the underlying point process satisfies the hypothesis of Lemma 15 


4.2.1. Definition of the good boxes. For K > 1 to be determined later, consider a partition of 
boxes of side K: 


into 


B z = B? := Kz + 


K Kid 

-,— . z E 7L d . 

2 ’ 2 J ’ 


We say that a box B Zl is nice if: 

-i- #(B zl nf) <c 2 K d , 

-ii- when B Zl is (regularly) cut into := (4[~\/cf| +5) d sub-boxes fr Z| of side s := K/otd, each of these 
sub-boxes contains at least one point of £. 

For B Zl and B Z2 two neighboring boxes, we write 6 Zl,Z2 , i = 1, ..., + 1, for the sub-boxes of side s 

intersecting the line segment [Kz\,Kz^[. We denote by c Z1 ’ 22 the center of the small box b Zl,Z2 and by 
u Zl,Z2 the nucleus of the Voronoi cell of c Zl,Z2 . 

We finally say that the box B Zl is good for G n if it is nice and: 

-Hi- for any Z 2 £ Z rf with ||zi — Z 2 11 00 = 1 and any i £ {1 ,..., 2[ \Vd] + 3}, M G c *i,* 2 ( 1+V Q) S holds. 


4.2.2. Construction of paths. For each good box B Zi , we choose the nucleus of the Voronoi cell of Kzi as 
reference vertex and denote it by v z% . Thanks to -ii- in the definition of good boxes, one can check that 
v Zi £ B Zi . 

Consider two neighboring good boxes B Zl and B Z2 . We want to verify that there exists a G n (£) _ 
path between v z , and v Z2 so that condition ([ 2 ]) of Proposition [4] holds. To this end, let us consider the 
collection u Zl,Z2 = v Zl ,.. ., u Z2,Z2 i = v Z2 of the nuclei of the Voronoi cells of the centers of the sub-boxes 
6 Z1 ’ Z2 , 1 < * < oid + 1- We claim that, for i = 1,..., cc^, there is a G n (£)-path yf 1 ’ 22 between u Zl,Z2 and 
u z !f 2 included in B Zl U B Z2 . Condition -ii- in the definition of nice boxes implies that u Z1 ’ 22 and u z _J _’ Z2 
belong to c Zl,Z2 + [—(1 + Vd)s, (1 + Vd)s] d . The existence of qf 1 ’ 22 then follows from -in- in the definition 
of good boxes and the fact that c Z1 ’ 22 + [—(1 + Vd)s, (1 + Vd)s\ d C B Zl U B Z2 . Finally, we obtain the path 
between v Z] and v Z2 needed in Proposition [4] by concatenating the y^s and removing any loop. This path 
has length at most 2 c 2 K d since it is simple and contained in the union of two good boxes. 


4.2.3. Domination. Since the underlying point process f is supposed m-dependent, the process of good 
boxes has a finite range of dependence k. Clearly, the fact that a box B Zl is nice or not depends only on 
the behavior of £ in this box. Moreover, due to the definition of the edge set of G n (£), the fact that two 
vertices u,v £ £ are connected by an edge in G n (£) only depends on B(u, n||« — u||) fl£. Hence, in order to 
decide if c *i ,*2 (i+Vd) s holds or not, we only need to know (c zl ’ Z2 , ( 2 n+ + Vd)s) Fl^. Thanks 

Theorem 0.0], we know that there exists p* = p*(d,k) < 1 such that if V[B Z is good] > p*, 


to 


LSS97 
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then the process of the good boxes X = {X z , z £ Z d } dominates an independent site percolation process 
on 7L d with parameter greater than p® lte (Z 2 ). We will chose s and thus K such that V\B Z is good] > p*. 
Using (V) and (D), one can show that, for s large enough: 

\ — v* 

V [ B z is not nice] < —-—. 


Moreover, using stationarity and isotropy of £ and Lemma 15 


V[3z ~ z ,3i £{ 1 ,... ,2\Vd\ + 3} s.t. A ( , ^y (\ + ^d)s ^ oes no ^ 

2|" v / d]+3 

^ S S r ^ A G a ,c^ 2 ,(i+vd) S y] 


z' ~z i= 1 


= M(2lVd\+3)V[(A aMl+v - d y] 

= 2d(2\\/c[\ + 3 )ci 5 exp ( - (1 + Vd) ci6 s cw ) 

1 -p* 


< 


for s large enough. We have just verified that the hypotheses listed in Subsection |1.2| allow us to apply 
Propositions [2] and [4] in the cases of Delaunay triangulations, Gabriel graphs and creek-crossing graphs. 
Theorem Q] thus follows. 

Appendix A. Technical results for the proof of Lemma |TT] 


In this appendix we collect preliminary results used in the proof of Lemma 11 and the assumptions of 
Proposition [4] are supposed to hold. 

A.l. Stationary point processes and their Palm versions. We first recall useful facts concerning 
stationary point processes and their Palm versions. 

Lemma 16. (1) For p >1, if A is the disjoint union of k unit cubes, then 

T [# (£ n A) p } < k?e [#(£ n [-1/2, l/2 ] d ) p ]. 

(2) Let 1 < p < 8. Let (A n ) ne pij be a convex averaging sequence of bounded Borel sets in that 
is a sequence of convex sets satisfying A n C A n+ \ and r(A n ) — > oo where r{A) := sup{r : 
A contains a ball of radius r}. The following ergodic result holds: 


#(£nA n ) a.s., LP(Af,V) 


AVol {An) 

where A denotes the intensity of the point process. 

(3) Let 1 < p < 8 . For any C £ M: 

AVo1(C2At) L p (At,V) 


> L 


#(£ n c 2 n) + CN d 1 


A 1. 


(4) For any h £ L 2 (Afo,Vo): 

£ [J 

(5) For any h £ L 2 (Mq,Vq): 


h(r x f)f(dx) = AVol(A)5 0 [/i], VA £ B(R d 


lim T[m^(/r(r x £))] = £ 0 [h], 

oo 


with xxv N as defined in Subsection 


3.3.1 
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(6) Fix 0 < a < 1 and define M = M(N , a) := N — |_iV"J. For any h £ L 2 (Mo,Vq): 


lim £ 

iV—>■ oo 


1 

#(£ n CW) + #(rjv) 



^-('TxOC(dx) 


= £o[h]. 


(7) //i/ie stationary point process satisfies (V) and £[#(£n[—1/2, l/2] rf ) 2 ] < oo, then its Palm version 
satisfies (V’)- 

(8) If the stationary point process satisfies (D) and £[#(£n[—1/2, l/2] d ) 2 ] < oo, then its Palm version 
satisfies (D’). 

Proof. Since the point process is stationary, 0 is immediate. Point ([2j is a consequence of [DVJ08| 
Proposition 12.2.II] and of the ergodicity of V. Observe that: 


£ 


\\o\{C 2N ) 


#(£n cw) + chvrf-i 


-1 


/ AfV \ p ^ 

s'-c)” 


+ £ 


#(£ n c 2 n ) < C7iv d 

AVol(CW) 




P 

1 1 #KnC 2J v)>C7A rd 


Thanks to assumption ([3]) in Proposition [4j the first summand vanishes in the limit. Using point ([2]), one 
has: 


/n(0 ■= 


AVol (C 2N ) 

#(£ n C 2 n) + CiY d_1 


p 

1 #(«nC' 2 jv)>c 7 Af d 


■ -> 0 
TV—>-oo 


Since |/jv| is bounded, the second summand goes to 0 by the dominated convergence theorem and point 
([3]) is proved. Point Q is a special case of Campbell formula (see for example SW08, Theorem 3.3.3]). 
Thanks to point Q, in order to derive point ([ 5 ]), it suffices to show that: 


"#(£ C C 2 n) + #(Tat) A Yol(C , 2 JV- 2 r c ) ' Jc 2N -2r c 


h(T x £)€(dx) 


= 0. 


lim £ 

AT—>-oo 































ANNEALED INVARIANCE PRINCIPLE FOR RANDOM WALKS ON RANDOM GRAPHS IN R d 


25 


Using four times the Cauchy-Schwarz inequality, and points [7| and[7}, one has: 


£ 


#(£ n C 2 n) + #(T N ) A Vol(C 2 Ar-2rJ ' Jc 2N - 


h{T x Q€( dx) 


2 r c 


<£ 


/ A Vol(C , 2 Ar- 2 rJ _ 2 #(■$ H CW- 2 rJ 

vn c 2 n) + ^(rw) / AVoifCh 


5 


A Vol(C , 2IV-2r c 


A Vol(C , 2AT_2r c ) 


#(£ n C 2N -2r c ) \Jc 2 n 


( / ^('TxOC(d^)) 

J Con— o.r„ 


<£ 


A Vol(C2jy-2r J 

#(£ n CW) + #(rjv) 


A Vol(C , 2AT_2r c ) 


5 


’C 2 N- 2 r c 


1 N | 2 #(C H C 2 N-2r c ) 

/ A Vol(C' 2 AT_2r c ) 
^(U:0 2 C(dx) 


=5 


s x £ 


£ x £ 


A Vol(C , 2AT-2rJ _ 

#(£ n CW) + #(rjv) 

/ A Vol(C , 2V-2r e ) 

'#(£ n C 2 n) + #(Tn) 

A Vol(C' 2 IV-2r c ) 


n C 2N -2r c ) 

/ A Vol(C' 2 IV-2r c ') 


£o[h 2 


- l 


-1 


£ 


#(£ n C 2N -2r c 


So [h 2 


Vol(C 2 7 V- 2 rJ 2 

£ [#(£ n [- 1 / 2 , l/ 2 ] rf ) 2 ] %o [/i 2 ]. 


'#(£ n cw) + #(IV) 

This completes the proof since h E L 2 {Mq,Vq), £[#(£ fl [—1/2, l/2] d ) 2 ] < oo and by point ([3j): 


5 


A Vol(C , 2 Af-2r e ) _ 1 N 4 
#(£n CW) + #(Tat) 


TV—s>oo 


A 0. 


The proof of point ([h]) is very similar to the previous one. 

Let us show that (V) and £[#(£ fl [—1/2, l/2] rf ) 2 ] < oo imply (V’). By the definition of the Palm 
measure, we have: 


Vo 


# (^nft + l-iAirf)) =o 


\ \ 

1 

) ) = 0 

= -£ 

)) 

A 


E 


L #(T^n(z+[-L/2,L/2] d ))=0 


yet n[-i/2,i/2] d 


where A denotes the intensity of the point process. Let us cut x + [—L/2, L/2] d into 2 d cubes of side L /2 
having x at a corner C l L / 2 , ■ ■ ■, C 2 • If the sum in the last expression is not 0 and L > 1, then at least 
one of the C^ 9 s is empty. Hence, the union bound and the Cauchy-Schwarz inequality imply that 


Vo 


# (e° n (x + [—L/2, L/2] d )) = o] < -£ # (e n [-1/2, l/2] d ) 

# (^n [- 1 / 2 , i/ 2 ] d 


#(^n(a:+c/y 2 )=0 for some je{1,...,2 d } 


d 

25 „ 
£ - £ 


p 


#(£n[-L/4,L/4] d ) = 0 
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Point 0 then follows from (V). Similarly, we can write that 


Vo 


#({°n(i + [-L/2,L/2]' 1 )) >c 5 L“ 


= a £ 

S / £ 


E 


L #(r^n(x+[-L/2,L/2] d ))>c 5 L‘ i 


_y&i n[-l/ 2 ,l/ 2 ]d 

# n [-1/2, 1/2]°^ l#(^ n (x+[-I,,L]<J))>c 5 L d 




# (?n [- 1 / 2 , i/ 2 ]' 

This proves point (M) thanks to (D) by choosing C 5 = 2 d C 2 - 


V 


#(£n[-L,L] d )>c 5 L a 


□ 


A.2. The mean boundary contributions are negligible. 


Lemma 17. Let h^ N : 


be a function verifying: 


|/4( x )| ^ c i7 


;,r c )) ifxeQ,^ 


N 

ifx£T N 


for some cyj < 00 independent of N. 

Then the mean boundary contributions vanish in the limit: 


lim £ 

N—±oo 


#(£ n c 2 n) + 


e 

X ^n\&N- T o 


|^(*)| ? 


= 0, 1<P<4. 


(42) 


(43) 


In particular, this is the case for h^ N = ip^ N or ip^ N defined in (28) 


Proof. Note that with (42): 


E 

X ^n\^N-t c 


But, on the one hand: 


14^)r < (P M* Tn J0^ + E n B2(x,r c )Y 


%-r c 


£ 


#r^ r) P 

< f 

#(1V) 2 

2 

£ 

r#(-Bh) 2p i 

#(? n C 2 Jv) + #T N #{£n) P 


_(#(£ n C’ 2 n) + #r tv) 2 . 



# r A 


N 


_(#(£ n CW) + #iv)‘ 


< Ci 8 5 


(44) 
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where cis is a constant which does not depend on N thanks to the point © of Lemma [T(i| On the other 
hand, using the Cauchy-Schwarz inequality and Lemma ©10: 


£ 


#(£ n c 2 n) + #r n 


#(^n J B 2 (x,r c )) J 


< £ 


#(£ 0 C 2 N \ C 2 N-2r c ) 

_ (#(£ n c 2N ) + #r N ) 2 

x £ 1 


#(£ O C 2 n \ C 2 N-2r c ) 

< (A Vol (C 2 N \ C 2 N-2r c )) 2 £ 
Vol {C 2 N \ C 2 N-2r c ) 


#(£nR 2 (x,r c )) J 


< 


Vol (C 2JV ) 


£ 


#(£ O C 2 n \ C 2N-2r c ) 

_(#(£ n c 2N ) + #rAr ) 2 

A Vol (C 2N ) 


£n 


#(£nR 2 (o,r c )) 


2p 


£n 


#(^B 2 (0,r c )) 


2p 


#(£ n c 2 n) + #r iv 

Hence, (43) is proved since the r.h.s. of ( |44[ ) and the r.h.s of (45) go to 0 thanks to Lemma 16 ([3]) . 
Let us check that ip^ N defined in |28| satisfies (42). If x E 0.)y, then: 

|Vat( x ’)|— l{a;,j/}eS G ( ? ),||x-j/||<r c |?/1 ~~ x l| 


(45) 




If x E Tat, then: 


|</4( x )| < 


+ #(r^) ^ + ^ + 

7ry J yerjv 

<(#(?0 B 2 {x,r c )) - l)r c + r c = #((nB 2 (i,r c ))r c . 

=# (-Bat) 


#(iv) 


X] |yi + -v| + 


2 ^ 


#( r v) 


X] k 1 ~ ^1 - 


y&B^+ 


#( r v) 


Hence, satisfies (42) with cyj = r c . The same computations show that ijfi N satisfies (42) with cn = r%. 

□ 

A.3. Control of exit times. Let us give a sketch of proof of the following result which is a slight 
rewriting of FSBS06, Lemma 5]. 

Lemma 18. Let and Tn,x be as in (33) and (35), 0 < a < 1 and M := N — |_IV“J. Then, 


lim £ 

N—>00 

lim £ 

N —>00 


m 


N 


(lxeC2M P Ar,3;[^7V — ^]) 


m %(^-xeC 2 M^ > r x ^[TN,x < ^]) 


= 0, 
= 0. 


(46) 

(47) 


Proof. Recall the notations from Subsections |3.2| and |3.3.1| We first note that the expectations in 
and (47) coincide. Actually, for each IV e N*, each locally finite set and each one can define 

a probability measure /a on x S° such that: 

KA x 3°) = P ^ X (A), VA e B(S 4), 

xB) = P Tx z(B), VB € B(3°), 
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and such that, /j-a.s., 7^ = Tn,x and ut = x + X t (£°), 0 < t < T^. The existence of such a coupling 
implies that P^xt^iv < t] = P t x ^[Tn,x < t] and we only need to prove (46). 

To this end, let us part C 2 N- 2 r c \ C 2 M into unit cubes Qi = Qn,i, i £ I = In- For n G N*, we set 
I* := {[:= (l\, - ■., l n ) G I n '■ lj A h, 3 7 ^ k}. For u G such that 7^ (w) < 00 , let us define: 

k = k(uj) : = #{i el: 30 < t < 7^(w), u t G Q*}, 

and i := ifc) G the list of such indices sorted in the chronological order of first visits. For 

j = 1,... ,k, we denote by Xj the first site of Qj visited by u and by tj the time spent at Xj during the 
first visit. 

Now let {if, i e I, £ e J\f} be a family of independent exponential random variables with respective 
parameters #(£ n Qi), where 

Qi := {y e R d : dist(y, Qi ) < r c }. 

Remark that given £, k and x = (aq, .. . ,xQ, tj (1 < j < k ) are independent exponential variables with 
respective parameters no larger than #(£n Qi). Thus, n, y being fixed: 


N,x 


[7^ < 11 k = n, x = y\ < + • • ■ + if < t]. 


(48) 


Since k > k m ; n := [N a \ — r c , one has using (48): 


£ 


m N Pat ,x [7/V — *]) 


#1 


= X X 5 [ m N^xeCiM X P NA 7 N- t > k = n '2 = y\) 

n=lN a ]-r c iei* yen^iQij^i 

#1 

< X] X 8 [ m 5r ( 1 ^eC 2M P^ ) . T [k = n,i = l]) x P [l* + ■ ■ • + < i] 

n=yN a \-r c l&I* 

It remains to bound the r.h.s of (|49|). To this end, for k > 0 and l G I*, let us define the event: 


(49) 


A = A(k, l) = {{eM: #{j: l<j< n, #(£ n Qj) > k£ [#(£ D Qi)] } > 

Due to the Markov inequality and the stationarity of V, one has: 

V(A) < —£ #{j : 1 <j<n, #(f n Qj) > k£[#(€ n Qi)]} 


n 

< 2V 


#(£n Qi) > k£[#(£ n Qi)\ 


a 0 . 


Note that if £ 0 A, at least |"n/2] of the exponential random variables have parameters 

smaller than n£[#(£ n Qi)]. It follows that on A c : 

(*£[#({ nQJjty 


P [T* + • • • + < t\ < p(n, k) := e ~^[#(^Qi)]t ^ 

r=Tfl 

Finally, we deduce by collecting bounds that: 


v\ 


£ 


m 


n{^ x ^ c 2M^*N,x \7n — *]) 

#1 

<£ 


^(ixeCaM- X ^ P N,x[ k = n ’i = -])( lA + lAcp (\- Na ^ ~ r ^ K )) 

n=[N a \-r c l_£l* 


< 2V 


#(^nQi) > «£[#(£ nQi)] + p{[N a \ -r c ,n). 
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This gives the result by letting N and then k go to the infinity. 


□ 


Appendix B. Examples of point processes 


In this appendix, we check that the assumptions given in Subsection 1.2 are satisfied by PPPs, MCPs 
and MHPs. 

B.l. Poisson point processes. Stationarity, isotropy, ergodicity, aperiodicity and the finite range of 
dependence condition directly follow from the definition of homogeneous Poisson point processes. It 
is also known that stationary Poisson point processes are almost surely in general position. Standard 
computations give (V), (D) and that ^(£n [0, l] rf ) admit polynomial moments of any order. The second 
factorial intensity measure of a PPP(A) admits A 2 as density w.r.t. the (2d)-Lebesgue measure (see e.g. 
|CSKM13[ (2.30)]). 

B.2. Matern cluster processes. Matern cluster processes (MCPs) are particular cases of Neyman-Scott 
Poisson processes (see CSKM13, p. 17B). Cluster processes are used as models for spatial phenomena, 
e.g. galaxy locations in space |KPBS + 99 or epicenters of micro-earthquake locations |VJ70|. 


MCPs are constructed as follows. One first chooses a PPP Y of intensity A called the parent process. 
For any y £ Y, a centered daughter process £ y is then chosen such that, given Y, {£ y } y cY are mutually 
independent PPP with intensity y in B(0,R). Then, £ := (j y&Y (y + £ y ) is a MCP with parameters 
A, y, R. It is clear that such processes are simple, stationary and isotropic. Since its parent process 
has a finite range of dependence and daughter processes have bounded supports, any MCP has a finite 
range of dependence. MCPs can be seen as Cox processes (see [CSKM13, p. 166]). Their (diffusive) 
random intensity measures y Y have densities Z 5 (x) := ^2 y£Y 1 B( y ,R)( x ) w.r.t. the Lebesgue measure. In 
particular, (d — 1) dimensional hyperplanes and spheres are y Y —null sets and MCPs are almost surely in 
general position and aperiodic. Assumptions (V) and (D) can be proved to hold as in |Roul5b, §6.2]. 


Standard computations also show that #(£n [0, l] rf ) admit polynomial moments of any order. The second 
factorial moment measure of a MCP(A, y, R) is given, for A\. A 2 C two Borel sets, by: 


A^(Ai xA 2 ) = £ 


E E 1 

sie^nAi X2G£nA2\{3:i} 


= £ 

= £ 


£ 


E E i 

xie^nAi x 2 e£nA 2 \{2i} 


Y 


/ Z l (xi)Z Y (x 2 )dxidx 2 
UA, Ja 2 

= I f £ [Z y ( xi)Z v (x 2 )] dxidx’ 2 , 

J J A2 

where we used the expression of the second factorial moment measure of a PPP driven by Z' 1 (see 
e.g. |CSKM13, (2.28)]) and the Fubini-Tonelli theorem. Hence, it is absolutely continuous w.r.t. the 
(2d)-Lebesgue measure and admits 


£ [Z Y ( Xl )Z Y (x 2 )} = y 2 ^A 2 Vol (B(0, R)) 2 + A Vol (B(xi, R) n B(x 2 , R))j 


as density. 


B.3. Matern hardcore processes. In 1960, Matern introduced several hardcore models for point pro¬ 
cesses. These processes are dependent thinnings of PPPs and spread more regularly in space than PPPs. 
Such models are useful when competition for resources exists (e.g. tree or city locations, see (Mat86| and 
references therein). 
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Given a realization £ of an initial PPP of intensity A and R > 0, let us define by: 

: ||® - 2/|| > R, Vy G £ \ {x}}. 

Then, £i is distributed according to a Matern I hardcore process (MHP I). 


One can show that this point process satisfies (V) in the same way as in Roul5b, §6.2]. Its isotropy 


and stationarity follow from the construction. Since the MHP I is stochastically dominated by the initial 
PPP all the other hypotheses can be directly deduced from the case of Poisson point processes. 

Given a realization £ of an initial PPP of intensity A, independent marks {T x } x ^ uniformly distributed 
in [0,1] and R > 0, let us define: 

£n := {x e £ : T x < T y ,My <E £n B(x,R)}. 

Then, £n is distributed according to a Matern II hardcore process (MHP II). One can check as above that 
MHP II satisfies the assumptions of Theorem [TJ 
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